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We consider some duality relations for models of non-interacting particles hopping on 
disordered one-dimensional chains. In particular, we discuss symmetries of bulk-driven bar- 
rier and trap models, and relations between boundary-driven and equilibrium models with 
related energy landscapes. We discuss the relationships between these duality relations and 
similar results for interacting many-body systems. 

§1. Introduction and general hopping model 

Among the simplest models of non-equilibrium statistical mechanics are one- 
dimensional transport models. For example, one may consider models of particles 
hopping on a chain, where a current is forced through a system either by coupling to 
reservoirs at the boundaries, or by forces acting in the bulk. Despite their simplicity, 
these models exhibit rich behaviour, and are the subject of ongoing studies.® 

Here, we are interested in duality relations between pairs of one-dimensional 
models. The properties of the related models may be quite different: for example, 
one may sometimes relate non-equilibrium models to equilibrium ones®~® or one may 
find mappings between models with different realisations of disordered rates.®'® - ® 
Our recent work has focussed on propagation of single (or non-interacting) particles 
in one dimension,®'® motivated originally by properties of glassy model systems.® 
However, such models have a broad range of applications, and have been the subject 
of many analytic studiesP® The purpose of this article is to comment on some rela- 
tions between the simple mappings that we have found and mappings in interacting 
(many-body) systems. In particular, a relation based on an inversion of the energy 
landscape occurs in several many-body systems®'® a well as in our analysis.® We 
aim to elucidate the origins of these mappings, particularly a duality between sites 
and bonds of Id chains: this is facilitated by studies of simple models for which the 
'heavy machinery' of many-body theory is not required. 

After reviewing our previous work, which concentrated on systems at equilib- 
rium, we discuss a set of disordered models where a uniform driving force acts in 
the bulk. We show how symmetries of these models result in a factorisation of their 
master operators that resemble the supersymmetric form used for time-reversible 
systems. We discuss the reasons for this, despite the breaking of time-reversibility 
by the driving force. Then, discuss duality relations between boundary-driven sys- 
tems and systems with conserved particles. These results are related to recent works 
by Tailleur, Kurchan and Lecomte.®'® 

We first define a disordered one-dimensional hopping model for non-interacting 
particles by specifying rates for hops from site i to sites i — 1 and i + 1, which we 
denote by £i and respectively. Let ni{t) be the density of particles on site i at 
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time t, with equations of motion 
d 

— rii(t) = £ i+1 n i+1 (t) + n-i7H-x(t) - {£{ + r^Jn^t). (1-1) 

for i = 1 . . . N. It remains to fix the boundary conditions. The simplest case is to 
use periodic boundaries, identifying site with site N. In that case, the equations 
conserve the total number of particles, YU=i n «W = ^> an d the equations of motion 
can be interpreted as a master operator as in Ref. [5]). However, we also consider 
an alternative case where we consider no and un+i as time-independent reservoir 
densities, allowing a boundary-driven system without a conserved density.® 

We will use an operator notation, exploiting the linearity of the equations of 
motion. We define a state \n(t)) = 5^i n «(£)K)> with a basis such that = dij for 
integer i. The equations of motion are m\n(t)) = W\n(t)) with 

N 

W = ^2\i){£ i+l (i + l\+r i ^ 1 (i-l\-(£ i + r i )(i\) (1-2) 
i=i 

(The interpretation of the states (0| and (N + 1| which appear in this operator 
depends on the boundary conditions, as described above.) 

§2. Duality relations in models with conserved density and periodic 

boundaries 

In this section, we review some previous results,-' restricting our analysis to 
periodic chains. In this case we have, writing W as to distinguish from the 

dual operator W^ 1 ^ 2 ^ below, 

N 

W W = [\i + 1) - |*» {n(i\ - + 1|). (2-1) 

2=1 

from which conservation of particles is apparent as ^2i{i\W^ = 0. This allows us 
to intepret as a master equation for stochastic motion of a single particle. 

This hopping model is dual to a second model which is of the same form, but with 
hopping on sites with half-integer indices. That is, one has densities ni (t),n,3 (t), . . . , 

2 2 

from which we construct a state \n(t)) = n i _i — 5) whose evolution is given 
by a master operator 

N 

wiw = EO i + !>-H-i» - Jl - L ^ + ID ■ ( 2 - 2 ) 

i=l 

Since the particles reside on half-integer sites, we associate integer indices with bar- 
riers between the sites. Thus, Li and Ri are rates for motion to left and right across 
the ith barrier. Then, is dual to if we take 
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Fig. 1. (Top) Illustration of an 'energy landscape', denned in terms of site energies Ei and transition 
state energies E i+ i. The rate for hopping from a site is related to the sum of the site energy 
(measured downwards as shown) and the adjacent transition state energy (measured upwards). 
(Bottom) On inversion of the potential, the transition state energies become energies of new 
sites with half-integer indices: transition rates between these sites are controlled by the energies 
Ei which now have an interpretation as transition state energies. From Ref. |8|). 



whose physical interpretation as an inversion of the energy landscape will be dis- 
cussed below (see also Fig. [I]) . 

To reveal the duality between the models, we write 



JV 

i=i 
N 

J = £l*+5>M i |-4+l< i + 1 0> ( 2 - 4 ) 



i=i 

so that 

W W = -DJ, W (1/2) = -(JD)l (2-5) 

This factorisation means that and VF^ 1 / 2 ) have the same eigenspectra. For 

example, if is a left eigenvector of W^ 1 ' with eigenvalue A ^ then D^\ip) is a 
right eigenvector of W^ 1 / 2 ' with the same eigenvalue. To express this duality in a 
more standard form, we could write 

W (1) D = -DJD = D(W i - 1/2) ) ] . (2-6) 

Despite the simple form of (|2-3|) . we emphasise that this duality relates distinct 
pairs of hopping models. To illustrate this, we parameterise the rates in a region of 
the chain as 

— Ei — E. i ~Ei — E. ,i , „ _> 

£i = e n = e l+ * (2-7) 
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where we interpret the E{ and E i , 1 as site and transition state energies. The duality 

condition (|2-3p can then be interpreted as a swap of transition state and site energies, 
or equivalently as an inversion of the energy landscape, as in Fig. [TJ Our sign 
convention is that the Ei are measured downwards from an arbitrary baseline, so that 
site and transition energies are dual to each other. We note that a parameterisation 
of the rates in the form (|2-7|) is always possible on any subsection of the chain, 
but applying it to the whole set of rates requires additionally a global constraint of 
detailed balance: Hi r i = Ili^- 

In addition, one may relate the propagators within the two models. Interpreting 
the equations of motion (|l-ip as a master equation for a single particle we identify the 

propagator G^ \t) = (i\e wWt \j) as the probability that the particle is on site i given 
that it was on site j a time t earlier. In the description in terms of non-interacting 
particles, the steady state two-point connected correlation function is simply 

{ni{t)rij (0)) ss - (ni) ss (nj) ss = Gij(t)(nj) ss (2-8) 

where we use the label 'ss' to indicate a steady state average. Starting from (|2-6|) . 
we consider the matrix elements (i\e w 1 t D\j — \) = (i\De^ w 1 \j — arriving 
at 

where G^f' = (i + ||e w(1/2 '*|j + |) is the propagator in the dual model. 

We emphasise that such results have application beyond single-particle models: for 
example, the same relation applies to models of diffusing and annihilating defects in 
random potentials® (with some restrictions on boundary conditions). 

Importantly, Equ. (|2-9p allows the propagator in the dual model to be calculated 
from that of the original model, without any knowledge of the disorder. If one 
then chooses the and ti from (different) distributions that are independent under 
translation in space then one may prove that the disorder-averaged propagators 
satisfy 

4 1} (*)=4 1/2) (*) ( 24 °) 

with both sides depending only on the difference i — j. 

Returning to models with detailed balance, = rii r i> we now define an 

operator E^> such that M 1 ' \i) = Ei\i) and similarly i^ 1 / 2 ) sucn that .E^ 1 / 2 ) | j _|_ I) = 
+ \)- Then, taking 

N 

S = -De-* (1/2) = J2 (\i + 1) - l*» (i + (2«H) 
i=i 

N 

S = tfe-* (1) = J2(\i + l)-\i-k)){i\e- Ei (^) 
i=i 

we arrive at a more symmetric factorisation of the master operatorsP 1 

= SS, = SS. (2-13) 
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We previously considered consequences of this symmetric factorisation,^ 1 and denned 
a renormalisation scheme that acts symmetrically on S and 5. For the purposes of 
this paper, the key point is that detailed balance ensures that wWe- E(1/2) = 
e -£ 1 2 ) (yy(i/2)^ ( anc i similarly for W^). This additional symmetry of the master 
operators allows the symmetric factorisation of (j2- 13j) . which also implies additional 
'duality' relations W^S = SW^™ and SW® = W^S. This structure appears 
in supersymmetric field theories P 1 >EJ)Il2l j n supersymmetric models the two symme- 
try relations are linked to time-translation invariance and time-reversal invariance^J 
(via detailed balance and the fluctuation-dissipation theorem.) 

§3. Bulk-driven pure trap and barrier models 

We now consider disordered systems with uniform driving forces applied through- 
out the system. This clearly breaks detailed balance and time-reversal symmetry. It 
is therefore somewhat surprising that we can identify a restricted set of such models 
for which a factorisation similar to (|2-13|) is possible. We will discuss how this fac- 
torisation reflects symmetries of the models under a change in the direction of the 
bias, at fixed disorder. 

In this section, we restrict the form of the disorder to 'pure trap' and 'pure 

barrier' models. The former are obtained by setting E-.i = for all i, but choosing 

&-t- 2 

the Ei freely. Then, describes a pure trap model with = £{ = e~ Ei : site i 

represents a trap from which particles hop to left and right with equal probability, 
but the overall rate depends on the trap depth Ei. We add a bias to this model by 
taking 

n = xe- Ei , £i = ye- E \ (3-1) 

We take x + y = 2 for convenience since this factor simply rescales time, and we 
assume x > 1 for concreteness. Pure barrier models are dual to pure trap models, 
so their sites have half-integer indices and their master operators are of the form 
W^ 1 / 2 ). In the absence of a bias, all sites have the same energy Ei = 0, but rates 
for hopping between sites depend on the barrier being crossed. In the presence of a 
bias, right- and left-going rates are multiplied by x and y respectively. 
The master operator for the trap model is constructed as in (|2-5|) : 

W rT = W {1) = -DJ xy (3-2) 

where the label 'rT' indicates that we consider a trap model biased to the right, and 
we indicate explicitly that the current operator J depends on x and y (as well as on 
the Ei). One can then construct a master operator, W m = —(J xy D)^ that is dual 
to W rT , in accordance with (|2-5[) . It has barrier-crossing rates Ri = li = ye~ Ei and 
Li = ri = xe~ Ei and is therefore biased to the left, justifying our notation W eB . 

As noted above, these biased models do not respect detailed balance, so they 
may not be factorised as in (|2T3p . However, let W iT be the master operator for a 
trap model with a bias to the left (that is, W iT = —DJ yx ). It may then be verified 
that 

W iT = S'S, W m = SS> (3-3) 
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Fig. 2. Illustration of (local) energy landscape for biased trap and barrier models, following Fig.[T] 
A landscape inversion as in Fig. [T] relates (a), a trap model biased to the right, to (b), a barrier 
model biased to the left. However, Equ. (|3-5[) means that (a) is also dual to (c), a trap model 
biased to the left. Finally, inverting the landscape (c) leads to (d), a barrier model biased to 
the right. Thus all four landscapes represent models with the same spectrum, as described in 
the text. We emphasise that the pairs (a,c) and (b,d) are related by an inversion of the bias, 
but with the disordered rates Ei remaining invariant. The energy E3 is identified for all four 
landscapes. 



where S was defined above, and 

N 

S' yx = ^2(y\i + l)-x\i))(i + ±\ (3-4) 
1=1 

Thus, W m has two alternative factorisations, as — (J xy Dy from (|2-5|) and as SS' yx . 
These relations show that it is dual both to W rT and W eT , and all three operators 
share the same eigenspectrum. By eliminating W^ B , one can then also write a simple 
duality between trap models with opposite biases, 

W tT e m =e m (W*y (3-5) 

This reflects a physical relation between time-reversed trajectories: a trajectory Tr(t) 
in the right-biased trap model and the time-reversed trajectory 7f(t) in the left-biased 
model have the same probabilities in the steady state. Note here that both models 
have the same steady state densities as an unbiased trap model, (nj) ss oc e Ei . 

In an analogous manner we can eliminate from the duality relations that 
link it to W rB and W iB to get a duality relation between left- and right-biased barrier 
models: 

W m SD = SD{W r y (3-6) 

where SD = D^e~ EW D. This equation also expresses a relation between propaga- 
tion in the models, under inversion both of the bias and of the direction of time. 
However, we did not find any simple physical interpretaton of this relation. 

To understand the relationships between the four biased trap and barrier models, 
it is helpful to think in terms of tilted energy landscapes as shown in Fig. [2J Of course 
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the bulk bias on an entire periodic chain cannot be represented in this way, but if 
we restrict attention to propagation over finite times t and therefore finite regions 
of the chain, this is not an issue. For example, to parameterise the rates (|3-ip of 
W rT in terms of site and barrier energies as in (j2-7|) . one may define E'- = Ej + jF 
and £". 1 = —(j + \)F — 5 where F = ln(x/y) is the driving force (in units where 

the temperature and lattice spacing are unity) and 5 = ^ln(x/y) is an offset to 
the transition state energies that simply acts to rescale the time in accordance with 
x + y = 2. The situation is illustrated in Fig. _%&). The factorization VF rT = —DJ xy 
leads to the dual model W m = — {J xy D)\ a pure barrier model biased to the left as 
shown in Fig. |2](b) . This model can alternatively be factorized as SS' yx , which has 
as its dual W eT = S' yx S, a trap model biased to the left (Fig. E|c)). 

These various symmetries also result in relations between propagators in the 
various models. If the propagator for one model (say the right-biased trap model) is 
known then those of its three dual models may be constructed. For example, (|2-9p 
reads in the notation of this section 

Gff(t) - Gg_ x (i) = Gf_ (t) - Gf_ hiH (t) (3-7) 

independently of both bias and disorder; the labels on the propagators are the same 
as those on the associated operators W. However, noting that W S' = S' xy W rB 
one may derive by a similar method 

xGQ +1 (t) - yGQ(t) = xGPij+i® - yGf + i_. + ,{t). (3-8) 

which is again independent of disorder but now depends on the bias. One may also 
use SW rT = W rB S to arrive at a relation that depends on disorder but not on the 
bias. However, this last relation reduces to (|3-7p if one notes simply that 

GfJ(t)e^ = G^(t)e* (3-9) 

which follows from (|3-5p . The latter relation between propagation in left- and right- 
biased trap models can also be written in a disorder-independent manner: for finite 
t and hence i — j, the propagators obey detailed balance with respect to the effective 
energy landscape defined above, i.e. G T ?e E 'i = G^e E i. This transforms our last 
relation into 

Gf(^=(f)^Gf(t) (3-10) 

which holds independently of the disorder as anticipated. 

Finally, using (|2-10p and taking the disorder average of (j3-8|) . it may readily be 
shown that the disorder-averaged diffusion fronts in the various models are related 
as 

cZ.it) = Gt 3 {t) = -Gf_.it) = S£(t) (3-11) 

Thus, in this section and the preceding one, we have shown how symmetries 
under time-reversal (detailed balance) or bias-inversion may be combined with the 
general relation (|2-5|) . leading to symmetric factorisations such as (|2-13|) and (|3-3|) . 
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These relations allow the propagators of the models to be related to one another. 
They also reveal an intuitively reasonable link between left- and right-biased trap 
models, and a somewhat more surprising relation between left- and right-biased 
barrier models. 



We now turn to the boundary-driven hopping models described in Sec. [TJ That 
is, the equations of motion are given by (|l-ip for i = 1 . . . N but we introduce two 
new time-independent reservoir densities, uq and rtjv+i- The model depends on 
the rates rp, ii, . . . , rjv, £n+i which may always be parameterised in terms of site 
and transition state energies as in (|2-7p . The dependence on the reservoir densities 
appears only through the combination nQrQ and njv+i^/v+i, or equivalently, through 
noe~ E ° and nj\[ + ie~ EN+1 . Since the particles are non-interacting, if we increase both 
no and njy+i by a multiplicative factor this simply increases the total particle number 
by the same factor. Further, if the two reservoirs are at different chemical potentials, 
noe~ E ° — riN + ie~ EN+1 ^ 0, there is a current proportional to this difference. 

Writing \n(t)} = n «(*)K) as above (but now with reservoir sites included), 

the equation of motion takes the form Jj|P(t)} = W d \P(t)) where W d is of the form 
given in (|l-2p . Since the sum in (|l-2j) runs only from 1 to N, we have (OlW 1 = 
(N + l|Ty d = 0, consistent with time-independent reservoir densities. It may be 
verified that the total density in the model is no longer conserved, so that W d is not 
a stochastic operator: X^(*l^ d ^ 0- 

4.1. Duality relations 

Recently, Tailleur, Kurchan and Lecomte^'® considered boundary-driven mod- 
els of interacting particles, and showed that they are dual to 'equilibrium' models that 
respect detailed balance. They further showed that this duality can be demonstrated 
for disordered non-interacting systems, in a many-body representation. Here, we use 
the methods of the previous sections to analyse the non-interacting case. Since we 
use a basis of one-particle densities instead of a many-body basis, our analysis is 
simpler. 

We begin by factorising W d as 



§4. Boundary-driven case 



W d = AS d S' 



(4-1) 



with 



N 




(4-2) 



N+l 




(4-3) 



i=0 



N 



A 



\i)(i\ = 1 - [0}(0| - \N + 1}{N + 1 



(4-4) 



i=i 
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In S , we identify site —1 with site N + 1: we imagine a periodic chain with the 
reservoirs sites being adjacent, but we forbid direct hopping between the reservoirs 
by formally setting e = 0. This makes the duality relations simpler since there 
are then equal numbers of sites i and transition states i + |. See Fig. [3l 



(a) Absorbing/Boundary-driven 




12 - N N + 1 




Fig. 3. Illustration of energy landscape for boundary-driven models and models with reflecting and 
absorbing boundaries. We show only regions near the ends of the chain, (a) The boundary-driven 
and absorbing models have the same disorder: the large circles on sites and N + 1 represent 
reservoirs for the driven case (operator W A ) and absorbing sites for the absorbing boundaries 
(operators W^). There is an extra transition state on site N + I to facilitate comparison of dual 
operators, (b) The model with reflecting boundaries (operator W 1 ). The landscape is obtained 
by inverting that for the driven/absorbing systems, following the procedure of Fig. [T] There is 
an isolated site N + §, with no hopping in or out. This ensures that all operators act in spaces 
of dimension N + 2. 

Then, we define two new operators 

W T = S d AS d , W & = S d S d A (4-5) 

These operators describe hopping for conserved particles, with energy landscapes 
illustrated in Fig. [3l They all have equal eigenspectra: if \ip) is a right eigenvector of 
W d with eigenvalue A ^ then S \ip) is a right eigenvector of W T and S d S is a 
right eigenvector of W 11 , all with the same eigenvector. These duality relations lead to 
the same conclusions as the many-body analysis of Ref. [3]) (for the non-interacting 
disordered case). Indeed, the many-body duality transformation considered there 
can be written compactly in terms of the matrix elements of S d and ^fSJ Whether 
any similar transformation may be exploited in the systems of interacting particles 
considered in Refs. [2]),[3]) remains a question for future work. 

4.2. Consequences of the duality symmetry 

The duality symmetry allows relations between propagators in the models to be 
determined. In particular, for the steady state in the boundary-driven model, we 
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have {ni(t) nj (0)) ss - (ni)*^)* = Gf j (t)(n i ) ss with G§(t) = {i\e wH \j). From the 
relation AW* = W d A, it follows that for 1 < i,j < N, 

G%(t) = G* J (t) (4-6) 

where Gfj(t) = (i\e WH \j) is the propagator for the model with absorbing states on 
sites and N + 1. That is, steady state correlations in the boundary-driven steady 
state of the model W d are equal to propagation probabilities between sites of a 
model with absorbing boundaries. Similar results for interacting models can also be 
derived!^* Physically, one notes that once a particle visits a reservoir, all correlations 
associated with it are lost, so the non-trivial correlations in the boundary-driven and 
absorbing models are the same. 

Further, the duality between models with reflecting and absorbing boundaries 
can be obtained by a factorisation, either as (|2-5p or as (|2-13p . In the notation of 
this section, (|2-9|) reads: 

= G%(t) - Gj^t) (4-7) 

where we used (|4-6p in the second equality and assumed 1 < i,j < N. Remarkably, 
this allows the two-point correlations of boundary-driven models to be constructed 
from the equilibrium correlations of a dual model (W v ) with an inverted energy 
landscape but no net current. 

§5. Conclusion 

We have discussed duality symmetries of non-interacting hopping models that 
are revealed by factorising their equations of motion. These factorisations indicate 
symmetries of the models. For example, if particles are conserved, master operators 
may be written as = -DJ, allowing a dual model W 1 - 1 ^ = -L>tjt to be 

constructed as in (|2-5p . If the models have additional symmetry, these may be 
revealed through alternative factorisations such as = SS for models obeying 

detailed balance, or W rT = S' xy S for driven trap models, leading to new dual models 
jy(V 2 ) = ss or w rB = SS' xy . These factorisations and duality relations allow 
relationships between propagators within the various models to be derived. 

For boundary-driven models without conserved currents, a factorisation of the 
master operator exists which reveals its duality with a model where particles are 
conserved and boundary conditions are reflecting, which we intepret as an 'equilib- 
rium' model. This then also allows a further relation to be derived between the 
boundary-driven model and one with absorbing boundaries. 

Several of these relations are known in models of interacting particlesP^ al- 
though derivations in such models require considerably more detailed calculations. 
Whether the methods presented here can be used to interpret and generalise such 
relations remains a topic for future study. 
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